REGULARITY FOR SOLUTIONS OF NON LOCAL 
PARABOLIC EQUATIONS 



HECTOR CHANG LARA AND GONZALO DAVILA 

Abstract. We study the regularity of solutions of parabolic fully noti- 
ce . linear nonlocal equations. We proof C a regularity in space and time 

and, under different assumptions on the kernels, C 1,Q in space for trans- 
lation invariant equations. The proofs rely on a weak parabolic ABP 
and the classic ideas of [§] and [TO] . Our results remain uniform as a — > 2 
allowing us to recover most of the regularity results found in [10] . 



1. Introduction 
In this paper we are interested in studying regularity for solutions of 

(1.1) ut-lu = 0, 

where / is a fully nonlinear nonlocal operator. The previous type of equa- 
tions arises naturally when studying evolution equations for purely discon- 
tinuous Levy processes. In this case / is given by the jump part (which is a 
linear operator, usually denoted by L), i.e. for u 6 C^°, 

(1.2) 

Lu(x, t) = P.V. / (u(x +y,t)- u(x, t) - Vu(x, t) ■ y\B 1 (y))d(i(y) 



= h™, / ( u ( x + V,t)- u(x, t) - Vu(x, t) ■ yxBt (y))dfi(y). 

The type measures we are considering are of the form d/i = K(x, t; y)dy for 
some kernel positive K even in the y variable, i.e. K(x,t;y) = K(x,t; —y). 
This assumption allows us to rewrite the operator without the principal 
value in the following way, 

Lu(x,t) = 5(u,x,t;y)K(x,t;y)dy, 

where 25(u, x, t; y) = u(x + y,t) + u(x — y,t) — 2u(x, t) is the second order 
difference in space of u at (x,t). 

More complicated equations appear in competitive games, in this case 
the nonlocal operator is given by and inf sup or sup inf combination of linear 

l 
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u t - sup inf ( c a ,p + ba.fi ■ Vu + / 8(u,x,y) "f n+1 dy ) = 0, 



operators, 

(1.3) (Inf-sup type) Iu(x,t) = inf sup L a ^u(x, t), 

P a 

(1.4) (Sup-inf type) Iu(x,t) = sup inf L a gu(x, t). 

p a 

A particular case of the previous ones are the extremal operators which we 
denote in the following way 

(1.5) (Maximal) A4^u(x,t) = sup Lu(x, t), 

Lec 

(1.6) (Minimal) M.^u(x,t) = inf Lu(x,t). 

In general, we are interested in operators satisfying a standard ellipticity 
condition, which allow us to control I(u + v) — Iu by extremal operators. 
We will give a precise definition in Section [2j 

In the case of a family of symmetric kernels K a p(y) = a a ^{y)\y\~^ n+l \ 
L. Silvestre studied in [8] the regularity of the solution of the Hamilton- 
Jacobi equation 

a afi(y) , 

p r'" J \v\ 

where c a ^ is a family of constants, 6 Qj( g is a bounded family of vectors and 
A < a a r < A. He was able to prove that the solution of the equation was 
classic using a non-variational approach to proof a diminish of oscillation 
lemma. A modification of the proof in [8] allows to get the regularity for for 
the same type of equations we study here, but the estimates would blow up 
as a — ^ 2. 

The variational problem was studied by L. Caffarelli, C. Chan and A. 
Vasseur in [2] by using De Giorgi's technique. Also recently, M. Felsinger 
and M. Kassmann in [6], obtained a Harnack inequality where the constants 
remain uniform as the order of the equation goes to the classical one by 
using Moser's technique. 

The focus of this paper is to study regularity of solutions of (II. 1|) , that 
remain uniform as a — > 2. This will provide a natural extension to part of 
the theory already developed by L. Wang in [TO] , 

The paper is divided as follows. In Section [2] we gave the proper defi- 
nitions of viscosity solutions and maximal operators. We will also specify 
what type of equations we are dealing with, write the hypothesis over the 
kernels and give some properties. For instance, we state there a comparison 
principle and the existence of solutions for the Dirichlet problem. Section [3] 
is the heart of this article. There we proof a very weak ABP type of esti- 
mate that will allow us to proof a point estimate in Section HI The strategy 
to prove the ABP consist in showing that we can cover the contact set of 
u (with respect to its convex envelope V) by a set where u does not sep- 
arate too much from T in a fixed portion. This will consist of two parts, 
controlling the convex envelope by above and below around a contact point 
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respectively. For the control by above, we iterate Lemma 5.1 in [8], for the 
second part we use an appropriated barrier to show that T does not decrease 
too fast. Section [4] deals with the point estimate which follows the ideas of 
|10| . Finally in Section [5] we state and proof our regularity results. 

2. Definitions and Preliminaries 

2.1. Non local operators. To be precise about the formulas we presented 
in the previous section, we need to ask and integrability condition to K 
around the origin, 

(2.1) [ \y\ 2 K(x,t;y)dy < oo. 

J Bi 

It allows us to write rigorously 

Lu(x, t) = J 5(u, x, t; y)K(x, t; y)dy 

not only when u, with compact support, is in C^°, but also when u(-,t) G 
C x ' x {x). 

Definition 2.1. u(-,t) G C ' (x), if there exists a vector v G W 1 and a 
number M > such that 

\u(y + x,t) — u(x, t) — v ■ y\ < M\y\ 2 for \y\ small enough. 

Notice that this definition implies \5(u,x,t;y)\ = 0(|y| 2 ) as \y\ is close 
to zero. This is why we can get rid of the principal value in the integral. 

We say that a family C of linear operators satisfy the integrability con- 
dition uniformly in fix [— T, 0] when the upper bounds in (12. 1|) can be taken 
independent of L G C and (x,t) G fi x [— T, 0]. 

We will consider as in [1] absolute continuous weights oo which measure 
the contributions of the tails to the non local operators. This allows to 
compute Lu even when u does not have compact support. 

Definition 2.2. The space of function L 1 ^) consist of all u : R n — » R such 
that 

IMIlih : = J \u(y)\u(y)dy < oo. 

Definition 2.3 (Non local operators). We say that I is a non local, fully 
non linear operator with respect to u, if for every u(-,t) G C 1 ' (x) D L {oj), 
Iu(x,t) is a well defined real number. 

For g G (0, 2) fixed, linear operators with kernels of the form K ~ 
|y|-(n+o-) or com bi na tioris of those (by taking supremums and infimums) are 
contained by this definition. In such cases we say that the operator have 
order a and use oj = 1/(1 + \y\ n+u ). 

We say that / is translation invariant in space if 
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where r is the shift operator 



T(x,t)u{y,s) = u(y + x,s 



+ t). 



2.2. Continuous operators. For time dependent problems, the natural 
topology to use in R n x R is the so called parabolic topology. It is generated 
by neighborhoods of the form B r {x) x (t,t — r], for a given point (x,t) G 
R n x R. For instance, in this topology a function / is continuous if and only 
if f(y, s) ->■ fix, t) as (y, s) ->■ (x, t~). 

Whenever we want to see that a linear operator L has "continuous co- 
efficients" we fix a smooth test function u in an open set O C R n x R and 
check if Lu evaluated at O is continuous. In the non local case, we need to 
use not only functions which are smooth in O but also that the contributions 
from their tails vary in a sufficiently smooth way. This is the motivation to 
introduce the following space. 

Definition 2.4. Let C(a, b; L 1 {oj)) be the space of function u : (a,b] — > R 
such that 

(1) for every t G (a, b], u(-,t) G L l {ui), 

(2) for every t 2 G (a, 6], ||tt(-,ii) - «(-,*2)||ii( w ) -> as h ->• ^ . 
/it comes additionally with the norm, 



The space of functions against which we test the continuity of / are given 
by parabolic second order polynomials and functions in C(a, 6; 

Definition 2.5 (Test functions). The space S = S(£l x (— T, 0]) of test 
functions is the set of all pairs (v,B r (x) x (t — r, t]) such that v G C(t — 
t, t; L 1 (a;)), B r (x) x (t—r, t]C!]x (— T, 0] and v restricted to B r (x) x (t—r, t] 
is a quadratic parabolic polynomial, i.e. 



Definition 2.6 (Continuous operators). We say that a non local operator 
I, with respect to ui, depends continuously on the position in Q, x (— T, 0] if 
for every (v, B r (x) X (t — T,t\) G S, we have that Iu is a continuous function 
in B r (x) X {t — r,t] (with respect to the parabolic topology). 

We can understand a little bit better how the space C(a,b; L 1 (w)) ap- 
pears as a requirement for the continuity of the operator in time. Without 
this condition, even the fractional laplacian would not be a continuous op- 
erator in B\ x (—1,0] with respect to any positive u. Take for example u 
equal to zero in B\ x (—1,0] and let vary u freely outside B\ x (—1,0]. 



V>\\c{a,b;V-{u))) = SU P \\u(-,t)\\ L l^y 



te(a,b] 



n 



n 




REGULARITY FOR SOLUTIONS OF NON LOCAL PARABOLIC EQUATIONS 5 



2.3. Ellipticity. In the classical stationary case ellipticity means that, for a 
solution u of a homogeneous problem, the positive eigenvalues of its Hessian 
control the negative ones and vice versa. Geometrically, the positive and 
negative curvatures of the graph of u control each other. A way to define 
this precisely is by imposing the following condition on F, 

M~(D 2 (u - v)) < F{D 2 u) - F{D 2 v) < M + {D 2 {u - «)), 

where we are using the notation from pp. 

When we have in mind equations of order a € (0, 2), we may want to use 
Co, the family of all linear operators L which are comparable to the fractional 
laplacian of order a (to be defined), in order to define the ellipticity of I as 
M.£ o (u — v) < Iu — Iv < M~£ (u — v). This family is however too big and 
in order to get further regularity we need to impose further assumptions on 
the linear operators. By this reason we give a definition of ellipticity which 
is more general. 

Definition 2.7. Let C be a class of linear integro differential operators. We 
say that a fully non linear operator I is elliptic with respect to the class C if 

(2.2) M c {u-v) < Iu- Iv < M^{u-v). 

Going back to the definition of Co = Co(A,a) (A > 1). The precise 
condition for L to be in Co with kernel K is the following one, 

( 2 - 3 ) (2-a)^—<K(y)<(2-a 



In this family the extremal operators take the explicit form 
sup ( 



A4^ Q v(x,t) : = sup (Lv)(x,t) 



_ r A5 + (v,x,t;y) - A M (v,x,t;y) 

J |y| n+cr 

M Co v(x,t) : = inf (Lv){x,t) 

f A~ 1 5 + (v,x,t;y) - A5~ (v,x,t;y) 

= ( 2 - J j^ps d v- 

i™ 

Where 5^ denote the positive and negative parts of 5 (5 = 5 + — 5~). 

The factors (2 — a) become important as a — > 2~ as they will allow us 
to recover second order differential operators. 

Holder regularity for the spatial gradient of u requires ellipticity with 
respect to a smaller class. Given po > we define C\ = C±(a, A, po) C 
Co(o~,A) by the family of operators L £ C\ with kernel K, such that 

/ •, f \K(y) -K(y-h)\ , 

(2.4) / ' yy> lhl W JA dy < A, 

Jm. n \B P0 \ n \ 
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for every \h\ < p /2. It is sufficient that \DK\ < A/(l + \y\ n+a )) for d23| to 
hold. An important property of this stronger condition is that the parabolic 
equations associated with them remain invariant under scaling. We will not 
bother to give this family a specific notation because we will not use it in 
this work. 

2.4. Viscosity solutions. For viscosity solutions of an equation ut — Iu = f 
we always assume a minimum requirement of continuity for u. Here we 
denote the space of upper semicontinuous functions in fl x [— T, 0], always 
with respect to the parabolic topology, by USC(Ct x [— T, 0]). Similarly, 
LSC{Cl x [— T, 0]) denotes the space of lower semicontinuous functions in 

n x [-t,o\. 

With respect to the time derivative, it is natural for the parabolic topol- 
ogy to consider only the values of u towards the past. In this sense 



Definition 2.8. A function u G USC(Q x [-T, 0]) (u G LSC(tt x [-T, 0]) ), 
is said to be a sub solution (super solution) to ut — Iu = f, and we write 
iH — Iu < f (ut — Iu > f ), if every time (v,B r (x) x (t — r,t\) G S touches 
u by above (below) at (x,t), i.e. 

(i) v(x, t) = u(x, t), 

(ii) v(y,s) > u(y,s) (<p(y,s) < u(y,s)) for every (y,s) G B r (x) X (t - 

T,t]\{(x,t)}, 

then v t - (x, t) — Iv(x, t) < f(x, t) (v t - (x, t) — Iv(x, t) > f(x, t) ). 

An equivalent definition holds if instead of using parabolic second order 
polynomials as test functions we use test functions (p with less regularity 
around the contact point. This is important when we want to prove the 
maximum principle by means of and inf and sup convolutions. We omit 
it here and just assume that the maximum principle for viscosity solutions 
holds. The ideas for the proof of this result are standard and can be found 
in p] or in the appendix of [8] for the non local case with a = 1. 

The following example illustrates the importance of having test functions 
in C(a, b; L 1 (w)). Consider u(x,t) = Xex{o} where E CC W 1 \ B\. In the 
domain B\ x (—1,0) the function u satisfies u t ~ + (— A)°"/ 2 u = in the 
classical sense. When t = the equation is not satisfied any more, u t - (x, 0) 
is still zero in B\ but (— A) a / 2 u(x, 0) becomes strictly positive in B\. If we 
consider now the same equation in the viscosity sense, u is a solution even 
when t = 0. The restriction for the test functions to be in C(— r, 0; 
(in the case the contact occurs at t = 0) implies that such test function will 
not be able to see that u has a jump at t = 0. 

2.5. Qualitative properties. Most of the qualitative behavior of solutions 
of fully non linear non local operators, as considered by us, have been already 
proven in [3] or in the appendix of [8]. Here we state some of the results 



u t - (x, t) = lim 



u(x, t) 



u(x, t — h) 



h 
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already known, that we will need to use later on. The first lemma was 
proven in [3] and says that for a super solution regularity by below implies 
that the operator can be evaluated in the classical way. The following three 
results are the expected maximum and comparison principles which can be 
proven as in [8]. The last theorem regards to the existence of viscosity 
solutions of the Dirichlet problem by Perron's method, here we show that 
by using and appropriated barrier the solution achieves the boundary values 
in a continuous way. 

Lemma 2.1. Let I be a elliptic operator with respect to Cq and f a contin- 
uous function. If we have a super solution, m — Iu = f in Cl x (— T, 0] and if 
is a C 2 function that touches u from below at a point (xq, to), then Iu(xQ,to) 
is defined in the classical sense and (p t -(xo,to) — Iu(xo) > f(xo,to). 

Theorem 2.2 (Equation for the difference of solutions). Let I be a contin- 
uous elliptic operator with respect to Cq and f and g continuous functions. 
Given u and v such that Ut — Iu < f and v± — Iv > g hold in Q x (— T, 0] 
in the viscosity sense, then (u — v)t — A4^(u — v ) < f — g also holds in 
ft x (— T, 0] in the viscosity sense. 

Theorem 2.3 (Maximum principle). Let u be a viscosity super solution of 
u t - M^ Q u >0mOx (-T, 0]. 

Then 

inf u = inf u. 

nx[-T,0] ((R"\n)x(-T,0])U(R"x{-T}) 

Corollary 2.4 (Comparison principle). Let I be a continuous elliptic op- 
erator with respect to C$, u be a viscosity sub solution and v be a viscosity 
super solution of 

w t -Iw = f inUx {-T,0}. 

Thenu< v in ((R n \tt) x (-T, 0])U(M n x {-T}) implies u < v m$7x(-T,0]. 

Existence and uniqueness of a solution in the viscosity sense follows 
from the comparison principle by using Perron's method. The additional 
ingredient we need is a barrier that guarantees that the boundary and initial 
values are attained in a continuous way. 

Lemma 2.5. Let a € (0,2). There exists a non negative function ip : 
R n x (-oo,0] -)• R such that: 

(1) V = in B x x {0} 

(2) i> t - -M+ (ct) V > m (R« \ Bt) x (-oo, 0], 

(3) V > 1 in (R n x (-oo,0]) \ (B 2 x [-k,0]), 
for some k universal. 

Proof. Let ip = <p(x) be the one from corollary 3.2 in [3] and 

k- 1 = inf \M + ip\ A 1. 

B a \Bi 
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Then ip = (ip — k 1 t) A 1 satisfy all the requirements. □ 

Theorem 2.6 (Existence). Let a E (0,2), Q be a smooth domain, I a con- 
tinuous elliptic operator with respect to Co and f and g bounded, continuous 
functions. The Dirichlet problem, 

u t - Iu = f inUx (-T,0], 

u = g in ((R rt \ O) x (-T,0]) U (M n x {— T}), 

has a unique viscosity solution u. 

Remark 2.7. The boundary data g only needs to be continuous at the points 
in ((R™ \ O) x (-T, 0]) U (W 1 x {— T}) with respect to the parabolic topology. 
With respect to smoothness of the domain, we only require that Q satisfies 
the exterior ball condition. 

Proof. Let u be the solution obtained by Perron's method, 
u(x, t) = inf {v(x, t) :v t - Lv > f in O x (-T, 0] 

v > g in ((IT \ Q) x (-T, 0]) U (M n x {—T})}. 

It can be shown that u £ C(Cl x [-T, 0]), solves Ut — Iu = f in Q x (— T, 0] 
in the viscosity sense and « > g in ((R n \ Q) x (-T, 0]) U (M n x {— T}), see 
[5]. We will show now that u attains the initial and boundary values by 
comparison with appropriated barriers. 

Let's see the case of initial values first. Let b : R n — > [0, 1] a smooth 
bump function such that supp(l — b) = B\ and 6(0) = 0. The function 
ijj(y,s) = b(y)+\\M^ b\\ oo s satisfies ip t -M^ Q ip > 0. Let (x, -T) G ttx{-T} 
and e > fixed. By the continuity of g, there exists 5 > such that 
\g(x,—T) — g(y,s)\ < e, given that \x — y\ + \T + s\ < 25. Consider the 
following barrier, 

fty, s) = g(x, —T) + e + 2\\gU ^ ^) + ^} . 

(3 it is constructed such that fit — M~c fi > 0. Let's see that /? > g in 
((M n \ n) x (-T,0]) U (n x {— T}) and° therefore g(x, -T) < u(x,-T) < 
/3(x, — T) = g(x, —T) +e according to the definition of u. If \x — y\ + \T+s\ < 
28, then j3{y, s) > g(x, —T) + e > g(y, s). If \x — y\ + \T + s\ > 25 then either 
\x — y\ > 5 and then b((y — x)/5) = 1 or (T + s) > 5 and then (T + s)/5 > 1, 
in any case 

( / y — X J* -\- S \ J* -\- s \ 

P(y,s) > -llslloo + 2\\g\\oo s tp ( — t — i J H — — \ ^ a(y,s). 

After having that u(x,—T) £ [^(x, — T), g(x, — T) + g] we use that e is arbi- 
trary to conclude that u(x, —T) = g(x, —T). 

Let's consider now the case of boundary values. Let (x, t) S dCl x (— T, 0] 
and e > fixed. Let 5 > such that |<?(x,£) — g(y,s)\ < e, given that 
\y — x\ < 5 and s € [t — n6,t], k is the one from Lemma 12.51 By making 5 
even smaller we can also assume that the ball B§u{x — (<5/4)ra) touches Q 
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by outside with n its normal vector. Let ip the function from Lemma 12.51 
The barrier, 

, * > . „ (y - (x - (<V4)n))4 
P(y, s) = g(x, t)+e + 2\\g\\ OB 1> ) 1 ' " ,s - t 



is constructed such that f3t—M~£ Q (3 > 0. It also remain above g in ((M n \fi) x 
(-T,i])U(Ox{-T}). If\y-x\ < 6 and s G [t-K^i] then/%,s) + 
e > 5(2/) s )- If (y 5 s) is outside the cylinder B$(x) x[t — k5, t], then it is also 
outside the cylinder Bg/ 2 (x — (5/4)n) x [t — kS, t], then -0( ^~( x ~W 4 ) n )) 4 ; s _ 
t) > 1 and /3(y,s) > — \\g\\oo + 21151100 > g(y,s). Then we conclude as before 
that u(x,t) = g(x, t). □ 

3. Partial ABP Estimate 

The classic ABP theorem says the following. If u satisfies ut — M~u > 
-f in B 1 x (-1,0] with u > in dB x x (-1,0] U B 1 x {-1} then, 



1 

n+l 



inf u~ <c\ (f + ) n+1 dxdt 
Bix (-1,0] \77{ u =r} / 

where the domain of integration {u = T} is the contact set of u with its 
parabolic convex envelope T. 

In the non local case there is no hope to obtain a similar result by 
integrating only over {u = T}. In fact, consider the function u(x,t) = 
({x^ — l)xs 2 ( x )- The contact set in this case has zero measure, however 
there is a constant C > such that ut + (— A) a u > — C holds in B\ x (—1,0]. 

To sort out this difficulty, we consider the set where u is between V and 
r + M, for some positive and universal M. The theorem we prove in this 
section is the following one. 

Theorem 3.1. Let f £ C([— 1,0]) positive and depending only on the time 
variable, po > such that 1/2 + 9y / n2 _1// ( 2_<J )po < 2 and let u satisfying 

u t ~ M C U ^ ~fXB 1/2 in B 2 x (-2, 0], 

u >0 in ((R n \Bi) x [-2,0]) U (R n x [-2,-1]}), 

sup u~ = 1. 

Bix(-1,0] 

T/ien, 

c < /I /(i) n+1 |{n(-, t) < r(-, t) + C^VM/^)} n £ 1/2+VH2 -i/(2-. Vo |^ 

/or some constants c and C depending only on n, A, o~o and po- 

Remark 3.2. TTie domain of the equation B 2 x (—2,0] can be reduced to 
B\ +£ x (—1, 0] and just indicates that we are going to need some room in the 
following proofs. 
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Remark 3.3. The radius po will be a fixed universal constant in future sec- 
tions. For this reason it may be noticed that we call also universal constants 
to some quantities that may also depend on pq. 

The following is a more general version of the same theorem. 

Theorem 3.4. Let SI C IR n a bounded domain and ^2 CC fl\ CC Oo CC f2. 
Let f G C([— 1,0]) positive and depending only on the time variable, d > 
and let u such that 



for some constant C , depending only on n, A, gq, f^j ^1, ^0 an d 

The idea of the proof is to cover {u = T} with a disjoint sequence of 
rectangles Kj = Qj x Lj, where Qj is a cube in space and Ij is an interval 
in time such that: 

(1) The measure of the union of the covering is at least a fix constant. 

(2) In a dilation Kj of Kj, u is most of the time between T and V + M. 

3.1. Preliminaries. Here we will fix the notation that we will carry on for 
the next results. 

From u we construct the following auxiliary functions which allow us to 
get important information for u. We enumerate them in the following list: 

(1) Let (xq, to) G B±X (— 1, 0] such that sup BlX (_ 10 ] u~ = \u(xq, to)\ = 1. 

(2) Let u(x,t) = inf s6 [_i >t ] u(x, s). 

(3) Let r(-,i) be the convex envelope of u(-,t) supported in B3. Notice 
that r is convex in space and non increasing in time. This is what 
we also call a parabolic convex function. 

(4) Let dT(x,t) C W 1 be the set of spatial subdifferential of T(-,t) at 
(x,t). Notice that for every t £ [-1,0], dT(B 3 ,t) = dT(B 1 ,t), there- 
fore every subdifferential p of T can be assume to be the slope of 
some supporting plane, namely x — > p • (x — xq) + h, to the graph of 
u(-,t). 

(5) Let h(-,t) : dT(Bi,t) ->Rbe the Legendre transform of r(-,t) cen- 
tered at Xq, 

h(p, t) = sup{/i : p • (x — xo) + h < T(x, t) for all x G B3}, 
or equivalently, 

h(p, t) = sup{/i :p • (x — xo) + h < u(x, t) for all x S Bi, 
p ■ (x — xq) + h < for all x G -B3}. 



u t ~ M c u ^ ~fXn 2 infix {-d a , 0], 

u > in ((R n \ n ) x (-tf 7 , 0]) U (M n x {-d a }). 



Then 




l/(n+l) 



sup 

f2 x(-c£ CT ,0] 



(6) 



Let $ = (dT(x,t),h(dr(x,t),t)). 
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This first lemma give some of the basic properties of the previously 
defined functions. 

Lemma 3.5. Let u, f , xq, T, h and $ as defined above and consider also 
At > 0. Then the following properties hold: 

(1) The domain of h(-,t) is non decreasing in time. i.e. dT(B\,t) C 
dT(B 1 ,t + At). 

(2) h is non increasing in time. 

(3) The function h restricted to {T = u} is Lipschitz in time. Specifi- 
cally, for (xi,t\) G {r = u} and pi G dT(xi,ti) 



Proof. The first two properties are consequences of the monotonicity of T. 
If at time t, the plane x — > p ■ (x — xq) + h is a supporting plane for the graph 
of r(-, t) then at time t + At it crosses or touches the graph of T(-,t + At) < 
r(-, t), therefore by lowering h we can find a supporting plane for r(-, t + At) 
with the same slope p. 



In order to see the next property notice that for every p G r( - ,i) and 

x G B 3 , 



This makes h(p,t + At) an admissible candidate in the definition of h(p,t), 
and therefore h(p, t) > h(p, t + At). 

For the second part, notice first that p\ G dT{B\,t\ + At) because of 
the first property, therefore Ah is well defined. We also have proved that 
Ah < in the second property, so assume that Ah < and consider the 
following test function, 



The infimum of u — v in W 1 x [—1, t± + At] is strictly negative and attained 
at some point (^2,^2) £ B\ x (ii,<i + At]. Indeed, the plane x — > p\ ■ (x — 
x o) + h(pi,t\) crosses the graph of u in B\ x {t\,t\ + At], otherwise Ah 
would not be strictly negative. 

We have that Vt(x2,t2) = Ah/2At and 8(v(-,t2),X2',y) > 0. This is 
immediate if X2 =t y lie both inside B\ or both lie outside B\ . If X2 + y G B\ 
we use that u~(-,t2) < 1 in order to see that \v(x2 + y, *2 ) I is at most 
2\v(x2, ^2 ) I - Because u~ < 1, the plane x — > p\ ■ (x — xq) + h(pi,ti) is above 
— 1 at some point in B\, lets recall also that the same plane is below zero 
in B3. It tells us that its slope is at most 1/2 and then v(x2,t2) — v{x2 + 
y,t 2 ) = -pi ■ y < \pi\\y\ < 1 < -^(^21*2)- Then, if x 2 - y is outside B\, 
S(v(-, ti),xi;y) = (v(x 1 + y, t{) - 2v(x 1 ,t 1 )) + v{x x - y, ti) > 0. 



Ah := h( Pl ,t 1 + At) - h{ Pl ,h) > -2||/|| ioo([tl)tl+At]) At. 



p • (x - x ) + h{p, t + At) < T(x, t + At) < F(x, t). 
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We conclude by the monotonicity of d/dt — MT^ , 

|i»([ti,ti+At]) < ("t - M~ c u)(x 2 ,t 2 ), 
< (v t - Mc a v)(x 2 ,t 2 ), 
Ah 

< 



2At 

□ 

3.2. Configurations of the covering pieces. In the following lemmas 
we will study how the solution u detach from T around the contact set. 
In order to keep the statements as simple as possible we will describe here 
some recurrent geometric configurations and fix the notation for them. For 
u satisfying the hypothesis of Theorem 13. 1\ k G N, At, po G (0, 1), (x\,ti) G 
B\ x (—1,0] we consider: 

(1) Rt = B n (xi) \ B ri+1 ( Xl ) for n = 2- l 2- 1 /(2-) /O0) 

(2) Si = Ri x [ti - At, ti - At/2}. 

Eventually (xi,ti) will be fixed to be in the contact set, k will also be fixed 
of the order of 1/(2 — a). 

Lemma 3.6. Let u satisfy the hypothesis of Theorem \ 3.1\ and for k G N, 
At G (0,1), (xi,ti) G B\ x (—1,0] let ri and Si as defined above. Given 
(pi, hi) G R n+1 , M > and ji G (0, 1) such that 

(1) T(x,t) >pi-(x- x ) + hi for (x,t) G B 3 x {h - At,h], 

(2) For every i = 0, 1, . . . , k — 1, 

\{U ~ (Pi • ~ XO) + hi) > ll/ + llL°°([ fl -AMi]) MT -i > n 5 »l ^ 

t/ien /or (x,t) G -B rfc+1 (x 1 ) x [t x — At/2,ti] we have 



u 



(pi • (a; - s ) + h) > II/ + ||i,»([ tl _ At)tl ])At 



if At G (0,r2) andMn(rl- a -r 2 k - a ) > K, for some constant K independent 
of a. 

Corollary 3.7. Let u satisfy the hypothesis of Theorem \S.l\ and for k G N, 

At G (0,1), (xi,ti) G Bi x (—1,0] let ri and Si as defined above. For every 
a < 2 there is some k ~ 1/(2 — a) such that if, 

(1) Ate (0,4), 

(2) (xi,h) G {T = u}, 

(3) pi G dT(xi,h), 

then there is some sufficiently small radius r G (0,ro), so that for S = 
(B r (xi)\B r / 2 (xi)) x (ti — At, ti — At/2] the following holds for every M > 0, 

\{u-(pi ■ (x - x ) + h(pi,h)) > ||/ + ||^([t 1 -At,f 1 ]) Mr2 } n ' g l K ^!i M -x 



\S\ 
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Proof of Lemma \3. 6[ Let /3 = (3q(x — xi/rf.) G [0,1] be a smooth bump 
function such that supp /3o = -B3/4 and Pq = 1 in Bi/ 2 - Moreover we can 
choose /3o such that MT^ /3 > if /3o(x) < j3\ for some positive constant /3\. 
We want to use a test function in B\ x [fi — At, ti] of the form 

= P(x) + m(t)/3(x) - ||/ + |U- ([tl _ At!tl]) (t - (t! - At)), 
P(x) = (pi • (x - x ) + /ii)xsi (^), 

such that m(t\ — At) = and m > 2||/ + ||^ C x>([ tl _At,t 1 ])At for i G [*i — 
At/2,h]. 

Assume by contradiction that, 

b !?£ At/2 , 1 ^- p ) < ll/ + H^a* 1 -A^i]) A *- 
-°'- fc+ i(^i) >< i ti_A */ 2 .*ii 

Then for some (x2,t2) G ^3r fe /4( x i) x [*i ~~ At,ti], 

inf (u — v) = (u — v)(x2,t2) < 0. 
R«x[ti-At/2,ti] 

Notice that i?3 rfc /4(xi) x [t± — At, t\] is contained in the domain of the equa- 
tion B2 x (— 2, 0] if At < 1 and ro < 1. We use Lemma 12. II in order to do the 
following computations on u at the contact point (x 2 ,t2)- We also use that 
5(u — v, X2,t2) = 5 + (u — v, X2,t2) and 5(P, x 2 ,t 2 ;y) > for every y E W 1 , as 
in the proof of Lemma 13.51 



- m'(t 2 )P(x 2 ) + m(t 2 )r- a M2 l3o(x2 - x ), 

< (ii t - M2 u)(x 2 ,t 2 ) - (v t - M2 v)(x 2 ,t 2 ), 

< (u - v) t (x2,t 2 ) - M2 (u - v)(x 2 ,t 2 ), 
<-A- 1 (2-a)[ 5+(u-vx2,t2;y) 

= -A-H2 -a)[ 6+(u- m p x 2 ,t 2 ;y) 

J\Ji3 \y\ n+a 

If x 2 + y G Ri, then |y| ~ r^, because X2 G -B 3rfe/ /4(xi). Also 5 + (u — m(3) > 
ll/ + IU° c ([i 1 -Ai,t 1 ])Afr? every time x 2 + y G Gj(t 2 ), where 

Gi(t) = {u - (p • (x - s ) + /ii) > ||/ + ||icx> ([tl _ At)tl]) Mr l 2 } n (ifc x {i}). 

Therefore we obtain, for some constant Co depending only on A -1 and the 
dimension, 



(3.2) 



k—1 

-m>(3 + mr^M Co (3 < -C (2 - a)\\f+\\ L ~ ([tl „ AtM) Mj2 i^^r?- 



i=0 
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As in [7], if m satisfies: 

m'(t) = d(2 - a)\\f + \\ Lao([tl _ At>tl]) Mj2 " <Wm(*), 

m(ti - At) = 0, 

with constants c\ sufficiently small and C 2 sufficiently large then we get a 
contradiction. Indeed, if fi§{x 2 — %i/rk) < A then .M£ A) > and we 
see that < c\ < Co implies a contradiction by substituting m' in (|3.2p . 
Otherwise, if /?o(^2 — xi/ru) > Pi, then we also get a contradiction in (|3.2p 
ifC 2 > \\M2 Po\\oo/Pi. 

We finally need to check that m > 2||/ + ||£oo([ il _Ai i t 1 ])At in [ti — At/2, t\] 
from the hypothesis of the lemma. We have an explicit formula for m, 

m(t) = Cl (2 - a)\\f^ {[tl ^ tl]) M^r^ f ^e'^^ds. 

Using the hypothesis (|3.ip of the lemma, for t >t± — At/2, we get 



m > C(2 - a)||/ + || Loo([il _ AMl]) M / iAt^__^_ e 



The quotient (2 - cr)/(l - 2 CT - 2 ) is bounded away from zero by a universal 
constant when a G [0,2]. Also e~ C2r k aM > e _C2 if At < r%. Finally 
in > 2||/ + ||L°°([ti-At,ti]) At is satisfied if Mii{r^ a — r^~ a ) > K for some K 
independent of a. □ 

The following is a geometric lemma that can be applied to any parabolic 
convex function. As a reminder, we say that T(x,t) is a parabolic convex 
function if it is convex in the variable x £ W 1 and non increasing in the 
variable t G M. 

Lemma 3.8. Let V : B3 x [— At, 0] — > M parabolic convex function such that 
\{T > M) n {B r \ B r/2 ) x [-At, - At/2] I 



\(B r \B r/2 ) x [-At, -At/2] I 



<£o- 



Then T < M in B r j 2 X [—At/2,0] if £0 is sufficiently small, depending only 
on n. 

Proof. By the convexity of T we can assume that its maximum N over 
B r/2 x [-At/2,0], is attained at (r/2e 1 ,-At/2). Therefore T > N in A = 
{(x,t) G (B r \B r/2 )x[-At,-At/2] :x-e\ >r/2}. Therefore, ife is smaller 
than \A\/\(B r \B r / 2 ) x [—At, — At/2]|, we obtain that N is necessarily smaller 
or equal than M in B r/2 x [-At/2,0]. □ 

By applying Lemma 13.51 and the previous lemma to T(x,t) — (pi ■ (x — 
x o) + h(p\,ti)) with all the hypothesis and conclusions of Corollary 13.71 we 
obtain the following result. 



REGULARITY FOR SOLUTIONS OF NON LOCAL PARABOLIC EQUATIONS 15 



Corollary 3.9 (Flatness of T). Let u satisfy the hypothesis of Theorem \3.1\ 

and (xi,ti), pi and r as in Corollary \3.7\ There exist a universal constant 
M > such that for every (x,t) G F = B r /2(x\) x [t\ — At/2,max{ti + 
At/2,0}] 

-2||/ + ||z/-([t 1 ,max{t 1 +At,0}])7' 2 S F(x,t) - (jp\ ■ (x - X ) + h(p U ti)), 

< M||/+|| L o 0([tl _ AMl]) r 2 . 

The previous corollary seems still insufficient to control 
|$(-B r/4 (a:i) x [t x - At/2, ti + At/2]) | 
\B r/4 (xi) x [ti - At/2, Fx + At/2] | 

The flatness property takes care of the n-dimensional size of dT{B r / i {x\) x 
{t}) by using the geometry of the convex function t) (t G \t\ — At/2, t\ + 
At/2]). Note also that the image of $(-,t) = (dT(-,t),h(dT(-,t),t)) is the 
graph of /i(-,t), for which we use again the properties in Lemma 13.51 

Corollary 3.10. Let u satisfy the hypothesis of Theorem \3.1\ and (xi,t±), 
pi and r as in Corollary \3. 1\ There exist a universal constant C > such 
that for K = B r / i {x\) x [t\ — At/2, max{ti + At/2, 0}] we have 

l*(*)l < C || /+Il n + 1 



|^| — Hi°°([ti-At,max{ti+At,0}]) - 

Proof. We do the proof for t\ + At < in order to avoid the difficulties that 
would arise if K goes beyond t = 0. In this case the proof does not differ to 
much from the one we present but makes the proof more technical. 

As a consequence of Corollary 13.91 we have that dY(B r u(x\) x {t\ + 
At/2}) C Bcr||/+|| i » ([tl _ Atitl] )(Pi) and then 

|ar(B r/ 40Ei)x{t 1 + Ay2}| + 

By Lemma 13.51 

§{B r/A (xi) x [ti - At/2,h + At/2]) C Cylinder, 

where 

Cylinder = {(p, h) :p G dV(B r/A {x 1 ) x {t\ + At/2}), 

h G [h(p,t),h(p,t) +2At||/+|| i co ([tl)tl+At]) ]}. 

Finally the measure of Cylinder is controlled by its base times the height. 

□ 

3.3. Covering of the contact set. We state now a weak version of the 
ABP estimate. The result consists in finding a covering of the contact set 
where the solution does not separate too much from the convex envelope in 
a given fraction of the union of the covering. 
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Lemma 3.11. Let u satisfy the hypothesis of Theorem \3.1l There exists a 
finite family of disjoint rectangles {Kj = Qj x Ij}, where Qj C W 1 is an 
open cube with diameter dj < ro/4 and Ij = (—(lj + l)At/2, —LAt/2) with 
lj non negative integer, such that: 

(1) Kj n {u = r} / 0, 

(2) \J j Kj^{u = Y}, 

(3) r is between two planes in Qj x Ij which are separated by a distance 

cwfh^d], 

(4) |^-)l<c|l/ll2iV,)l^^ 

(5) \{u < T+C\\f\\ LO o {h) d]}C\Kj\ > (l-e )\kj\, where Kj = W^EQj x 
[-(lj +3) At/2, -lj At/2]. 

Proof. Fix a slice B\ x I[ and cover it by a tiling of the form {Q x I{\ 
where Q have diameter ro/4. Discard all of those rectangles that do not 
intersect {u = T}. Whenever Q x I\ does not satisfy ([3]), (JI|) or ©, we 
split Q into 2™ cubes Q' of half diameter and discard all of the rectangles 
Q' x Ii whose closure does not intersect {u = T}. We need to prove that 
eventually all rectangles satisfy ([3]), ((4|) and ([5]) and therefore the process 
finishes after a finite number of steps. In fact we will show that it will finish 
before k ~ 1/(2 — a) iterations. 

As before, in order to avoid technical difficulties, we assume that I > 1. 
Let Qi X Ii D Q% X Ii D ■ • ■ 2 Qk X h 9 (#1, *i) such that (xi, t±) £ {u = T} 
and let's see that at least one of those rectangles satisfy all the properties 
([3]), ((!]) and ([5]). From Lemmas 13.91 and l3.10[ there is some radius r € [rfc,ro] 
and some subdifferential p\ € dT(x\,t\) such that the following are true: 

(1) |r - (p! • (x - xo) + /»(pi,*l))| < C\\f\\ Laa(ll) r 2 in F = B r/4 x [t x - 
At/2,ti + At/2], 

(2) MK)\ < CWfWl+^K] for K = B r/4 x fa- At/2, t x + At/2] , 

(3) |{u - (p! • (x - x ) + %i,ii)) > ||/|| Loo(/;) Mr 2 } n S\ < e \S\ for 
S = B r x [t x - At,h - At/2}. 

There is one of the rectangles Qj x with diam(Qj) = d, such that 
r/8 < d < r/4. Therefore Qj x I t C if(xi,ti) and conditions ([3]) and (jU) 
from the lemma are verified. To check ([5]), notice that S 1 C 16y/nQj x 
[— (Z + 3)At/2, —/At/2] = Kj, and that the volumes of S, Kj and Kj are 
comparable, hence 

\{u<T+\\f\\ Laa(Il) Md]}nK 3 \, 

> \{u <pi-(y- Xl ) + %i,ti) + \\f\\ L ~ {h) Md)} n S\, 
>{l-eo)\S\. 

This is how [i is chosen and this concludes the proof. □ 
3.4. Proof of Theorem 13. 1L We have as in [9] that 

$({r = u}) D Cone = {{p,h) : h G [— 1, 0], > 4|p|}. 
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The inclusion follows because for every (p, h) £ Cone the plane x — > p ■ (x — 
xq) + h can be brought from t = — 1 towards the future until it hits the 
graph of u (and also the graphs of u and V) for the first time. 
Therefore for some universal constants, 



c< i$(u^-)i < E i^(^)i <cE 

We group now the previous sum in each interval J\ = (—(/+!) At/2, —I At/2), 



i 



n+l 



£i* 3 -i<£ 



n+l 

L°°(Jl) 



U «i 



By Besicovitch, we can take a sub set of {Kj}i j= j t (denoted by the same) 
with the finite intersection property and still covering \Ji j= j t Kj such that 

E l^-|<cE l{«<r + C4- 1 /(2-)|| / || ioo{Ji)}n ^.|, 
ij=ji ij=Ji 

< c \{u < r + ca- 1 '^ \\f\\ L ~ {Jl) } n (uj^k, 



Notice that the contact set {u = T} can only occur where fxB P0 is positive. 
This implies that ( Uj j= j t Kj J C B 9 ^ po and then we have the following 



Riemann sum which is now independent of the covering, 



c < 



n+l 
L°°{Jl) 



{u<r + C4- 1 /(2-.)|| / || ioo(j;)} n B x j t 



where J/ = (—(I + 3)At/2, —/At/2). Now we just have to send At to zero 
to conclude the theorem. 

In [3] the partial ABP involves a Riemann sum of \f\ n which gets refined 
in the limit, when a goes to 2, and allows to recover the classic ABP. The 
estimate presented here is weaker, assuming that the right hand side of the 
equation is f(x, t) we notice that in Lemma 13.51 we need to take a global L°° 
norm in space and not just around the contact point. All the other proofs 
work fine in this sense. Our proofs can recover the following consequence of 
the classical ABP as a goes to two, 



l/n+l 



sup u < C 

BlX(-l,0] 



/ sup(/- 
'{r=-u} yeBi 



\n+l 



(It 



4. Point Estimate 

We are interested now in proving a point estimate that will allow us to 
control the oscillation of the solution. Our goal is the following theorem. 
(Recall that / G C[— 1,0] and non negative). 
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Theorem 4.1 (Point Estimate). Let <jq G (0,2) and a G ((To, 2). Suppose u 
satisfies 

U t~ M C (a) U ^ -/(*) in 5 1 X (- 1 ' ]' 

■1,0]. 



Then, for every s > 0, 
|{ U >s}n J B 1/2 x [-1, -1/2]| 



•e 



l-Bi/2 x [-1,-1/2] I V B i/2x [-1/2,0] 

for some constants e, C depending only on n, A and cq. 

The proof of Theorem 14.11 is done by induction as in [10J. The idea is 
to get a control of the measure of the set where u is bigger than a universal 
constant and then being able to reproduce the estimate at every scale. 

4.1. Initial configurations. 

Lemma 4.2 (Special Function). Let <to G (0,2), a G (o"o,2). There is a 
function p(x,t) G C(W l x [0,80]) and a constant C > 1, such that for any 
a G ((T , 2), 

Pt- ~ M c (a)P < - 1 + CXb 1/4 x(o,i] in B 4VE x (0,80], 

p(x, t)<0m ((R n \ B 2V7i ) x (0, 80]) U (R n x {0}), 
p(x,t) > 2 in Q 3 x [1,80]. 

The function p will also be C 1 ' 1 in the space variable and C 1 in the time 
variable (with respect to the parabolic topology), so that the computation of 
the equation is done in the classical sense. 



Proof. Consider 



inM n \5 a , 
in B s , 



where q is a quadratic polynomial chosen so that / is C ' across dB$. From 
Section 9 in [3] we know that it satisfies <M.£/ (T ^f > in W l \ B\j± for some 
sufficiently large p > and some sufficiently small 5 G (0,1/4), indepen- 
dently of a G (<7o, 2). By multiplying / by a sufficiently large constant we 
can also assume that 

M Co{a) f > 1 - C XB 1/4 in B 4y/E , 
\Df{x) ■ x\ < d in R n , 
f <C 2 in R n . 
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Consider also a continuous function m(t) > such that for some r S 
(0, 1) to be fixed, 

m(t) = t 1 ' 2 in [0,t], 

m(t) = t-VV 2 ^^ in (r,80], 

Let p(x, t) = m(t)f(y), for y = t~ 2 /°~°x. Notice that p(-, 0) = defines p 
continuously up to time zero. 

In the region B 4 ^ \ B±/ 4 x (0, r] we have that, 

(p t - M~ Co p){x, t) = ^~ {2/a )r 1 ' 2 Df{C 2 y) ■ y - t 1 ^ 2 ^ M~ c J(y), 

< {C 2 /2 + 2C l /a Q )t- 1 / 2 -t~ 3 / 2 , 

Therefore, by choosing r < l/(C2/2 + 2C\/oo) we make pj — .M^ p < in 
B 4VE \B 1/4 x [0,r]. 

In the region B 4y ^ x (r, 80] we have that, 

(p t - M Co p)(x,t) = m'(t)f(y) - ^Df(C 2 y) • y - ^M~ c J(y), 

<m'{t)C 2 + m{t)(C l + C Q )/T : 

which is zero by the construction of m in (r, 80] . 

Finally, we define p = A(p — B) + with B > chosen such that p — B < 
in ((M n \ B 2y ^) x (0, 80]) U (W 1 x {0}) and A > 1 chosen such that p > 2 in 
Q 3 x [1, 80] and p t - -M^ (ff) p < -1 in £ 4v ^ x (0, 80]. □ 

We are in shape now to prove a first control of the distribution. 

Lemma 4.3 (Base configuration). Let o~o £ (0,2), o~ G (o~o,2), an<i u a 
function such that 

u t -M~ oia) u > -1 in B^ x (0,80], 
u > in R n x [0,80], 

inf u < 1, 

Q 3 x[l,80] 

T/ien 

|{u > M } n Qi x [0, 1]| < /io|Qi x [0, 1]|, 
/or some universal constants fiQ € (0, 1) and Mo > 1. 

Proof. Let v = u — p, where p was constructed in the previous lemma. The 
graph of v goes below —1 at some point in Q 3 x [1,80], stays non negative 
in ((W 1 \ Q 3 ) x (1, 80]) U (R n x {0}) and satisfies 

v t ~ M C V - ( n * - M Co u ) ~ (P* - M C Q P) ^ -C'XBiMxfO,!]- 
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We apply now a Theorem 13.41 to v with f2 = B 4 ^, £Iq = Q3, O2 = B 1 / 4 
and Q,\ = Qi, 

c< f \{v < T v + C}C\B 1/2 \dt < \{u <p + C}C\Qi x [0, 1] I . 
Jo 

We just choose [1q = 1 — c and Mq = supg lX [ 0> i] P + C to conclude. □ 

We will need the following corollaries of the previous lemma. The first 
one will be necessary for the particular dyadic decomposition we will intro- 
duce in the next section. The second one iterates m times Corollary 14.41 



Corollary 4.4 (Flexible configuration). Let cjq G (0,2), a G (<to,2), r G 
[1, 8] and u a function such that 

U * ~ M C (a) U ^ - 1 m B ^ X (°> ( 3<T + M 

u > inR n x [0, (3 CT + l)r], 
inf u < 1. 

Q 3 x[r,(3- + l)r] 

T/ien /or /xi = , 

\{u> M }HQ 1 x [0,r]| < Ml |Qi x [0,r]|, 
/or /Uq € (0, 1) and Mq > 1 as in Lemma \4~. 



Proof. Notice that Q3 x [r, (3°" + l)r] C Q 3 x [1, 80] therefore if u goes below 
1 in Q 3 x [r, (3°" + l)r] then \{u > M } n Qi x [0, 1]| < /x |Qi x [0, 1]| and 

\{u > M } n Qi x [0,r)\ < /xolQi x [0, 1]| + |{u > M } D Qi x [1, r]| , 

<^|QiX[0,r]|. 

□ 

Corollary 4.5 (Iteration). lei cr G (0,2), a G (o"o,2), r G [1,8], > 1 a 
natural number, di = 3 ^zf^ an d u a function such that 

u > in R n x [0, d k r) , 
inf u < 1. 



Then 



>M k }r\Q 1 x [0,r]| < x [Ml, 

/or /ii G (0, 1) and Mq > 1 as in Corollary \4-4\ 

Proof. Just apply Corollary 14. 4| rescaled, k times. □ 
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4.2. A Calderon - Zygmund Lemma. The purpose of a Calderon - Zyg- 
mund type Lemma is to find a cover of a given set A with dyadic boxes that 
capture a fraction of A around a given /ii G (0, 1). 

The boxes are chosen from a dyadic decomposition of Q\ x [0, 1] that 
almost preserve the scaling of the equation. When a is either 1, log 2 3 or 
2 the decomposition can be made preserving the scaling of the equation by 
dividing by 2, 3 or 4 in time respectively. The following algorithm for general 
a was communicated to us by Luis Caffarelli. 

Initially we consider Q\ x [0,1], split Q\ into 2 n congruent cubes and 
take the 2 n possible cartesian products with [0, 1] to form the new dyadic 
boxes. In each step we consider one of the cubes Q r (xo) x [to, to + T a r\ 
and always divide Q r {xo) in 2 n congruent cubes; with respect to the time 
interval we do the following: 

(1) If r < 2 then we do not subdivide [to, to + r a r]. 

(2) If 2 < r < 4 then we subdivide [to, + r a T] in 2 congruent intervals. 

(3) If 4 < r then we subdivide [to, to + r a r] in 4 congruent intervals. 

Finally we take the cartesian product to form the new generation of dyadic 
boxes from Q r (xo) x [to, to + r a r\. 

This procedure verifies that if r G [1,8], then the boxes that Q r {xo) x 
[tojto + r^r] generates have dimensions r/2 (in space) and (r/2) <J r / (in time) 
for some r' G [1,8]. To prove it, just consider each of the cases. 

Given two dyadic boxes K and K we say that K is a predecessor of K 
if K is one of the boxes obtained from the decomposition of K. 

The following lemma follows as the one in chapter 4 of [I] with the 
difference that the Lebesgue decomposition theorem is applied to rectangles 
of dimensions p and p a instead of the standard cubes. 

Lemma 4.6. Let A C Qi x [0, 1] and p± G (0, 1), such that \A\ < p,\\Qi x 
[0, 1]|. Then there exists a set of disjoint dyadic boxes {Kj} such that: 

(1) | Uj Kj \A\=0, 

(2) \Ar\Kj\ > /J,i\Kj\, 

(3) \Ar\Kj\ < iii\Kj\. 

Proof. Starting with Q\ x [0,1], we subdivide the dyadic boxes (with the 
previous algorithm) that capture a fraction of A smaller or equal to /xi and 
select those boxes {Kj} that capture a fraction bigger than p\. Initially 
Qi x [0, 1] captures a fraction of A smaller or equal to fi\, therefore we know 
that Qi x [0, 1] is subdivided and Kj C Q 1 x [0, 1]. 

This process selects a family of disjoint boxes {Kj} that satisfy 2 and 3. 

To verify 1 we use the Lebesgue differentiation theorem. For each (x, t) G 

UjKj \ A there exist a family of dyadic boxes {K^ x '^ = Q r .(xj) x [U,ti + 
r f T i\}i>i such that, 

(1) (M)Gi^, 

(2) rj — > as i — > oo and Tj G [0, 8], 
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(3) lAnK^lKfulK^]. 

From K^'^ = Q ri (xj) x [ti,ti + r°"r] we construct a box with a scale a, 

Ki X,t) = Q Pl {xj) x [U,U + p%\ D K\ x,t) such that Pi = nr] 1 * . They satisfy 
instead, 

(1) (x^eR^, 

(2) pi — > as i — > oo, 

(3) \AnK^\ < MK^l with fix = i8 °°-^° < 1. 

Then we can apply a modified version of the Lebesgue differentiation theo- 
rem to conclude that | Uj Kj\A\ = 0. See for instance Exercise 3 in Chapter 

7 of nu. □ 

This lemma however can not be applied in our situation directly. The 
results of the previous section say that if u goes below 1 in some region in the 
future then we can control the distribution in the past; but the predecessor 
Kj might give no information of what happens with u in the future. For 
this reason we need to consider shifts in time of Kj. For a given cube K = 
Qx [to,to+r] and anatural number m > 1 let K m = Qx [to+r, to + (m+l)r]. 

The following lemma is proven as in section 3 of |10j. 

Lemma 4.7. Let A C Qi x [0, 1] and p\ G (0, 1), such that \A\ < pi\Q± x 
[0, 1]|. Then there exists a set of disjoint dyadic boxes {Kj} such that: 

(1) | U, Kj \A\ = 0, 

(2) \AnKj\ > m\Kj\, 

(3) \A\<^±^\Uj(Kj) m \. 

Proof. Select the same covering {Kj} of A from Lemma 14.61 Consider a 
disjoint sub covering of {Kj} which also covers A in measure (denoted by 
the same). Then 

\a\ < \A n Kj\ < mi Uj Kj\ < mi Uj {Kj u {kj) m )\. 
j 

Now we show that 
(4.1) | Uj (Kj U (Kj)™)] < ^±1| Uj (Kj) m \. 

Let E be the interior of Uj(Kj) m and consider the open bounded sets of the 
real line E x = {t £ 1 : (x,t) £ E}. Take the decomposition of E x into a 
countable set of disjoint open intervals E x = Ujl x and for every P x = (a, b) 
take TI X = (a (b-a), b). Define also TE X = UjTli and TE = U X TE X . 
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By Fubini, 

\TE\ = 



,i i - 



< I \E X \, 



m 

m + 1 , 
m 

m + 1 . . 
\E\ 



rn 



771 

On the other hand, the interior of L)j(Kj U (Kj) m ) is contained in TE, by 
the definition of T and the definition of the stacks. But L)j(Kj U (Kj) m ) 
and its interior have the same measure and therefore we obtain (|4.ip and 
conclude the proof. □ 

By combining Corollary 14.51 with Lemma l4.7l we get the following result. 

Lemma 4.8. Let o"o £ (0,2), a 6 (do, 2), m > itlq and u a function such 
that 

u t - M~ Co ^ a) u > -1 in B 2 ^ 3m x (0, d m ], 
u > in R n x [0, dm] 

Then for any i > 1, the dyadic covering {Kj} of A = {u > M[ +1 } n (Q\ x 
[0, 1]) with respect to the fraction fi± satisfies 

^(Kj) m C {u > M{} n (Qi x [0,d m ]), 

/or mo universal, M\ = M™, and \i\ G (0, 1), Mo > 1, as wi Corollary 

m 



Proof. Notice first that Mi = M™ implies by Corollary 14.51 that 

\A\ < \{u > M m } n Qi x [0, 1]| < ni\Q x x [0, 1]|, 

which allows us to apply Lemma 14.71 to A. 

By the construction of (Kj) m we know that L)j(Kj) m C Qi x [0, d m ] for 
m large enough. So we can assume by contradiction that there exists some 
Kj = Q r (xo) x [to, to + r a T] such that 

(1) vniffr yn u ^ 1> 

(2) i^nkji > /ui|ifj|. 

The rescaling (x,t) — >■ (r~ 1 (x — xo), r~ CT (i — to)) sends Kj to Qi x [0,r]; 
the stack (Kj) m to a subset of U^C^i x [diT, di + \r] if m is large enough 
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(at this moment we can fix mo); and it sends the domain of the equation 
Bn v /^ 3 m x (0, d m T] to a even bigger domain. 

The function v(x,t) = u ( rx+x ^ *+*o) sa ^ s g es 

V t ~ M C V - ~ - _1 hl X (°> dmT ^ 

v > in E n x (0,d m r], 
inf v < 1. 

U fcl Q3 a iX[diT,d i+ iT\ 

By Corollary 14.51 we get, 

\AHK 3 \ \{v>M^}C\Q l x [0, 1]| 

rp> ^r~rTi — Mi- 



|i^| " |Qi x [0,1]| 

But this contradicts the construction of Kj. □ 

4.3. Proof of Theorem 14. ll By combining the previous results we get the 
following discrete version of Theorem 14.11 which implies the proof of Theorem 
14.11 by rescaling and a covering argument. 

Lemma 4.9. Let <jq G (0,2), a G (o~o,2) and u a function such that 

u t - M~ C( u > -1 in B 2y ^ 3m x (0, d m ], 

u> inR n x [0,d m ], 

inf it < 1, 
Qix[C -l,C ] 

then for any k G N 

|{ U > Af 2 fc }nQ! x [0,1] I <^|Qi x [0,1] I, 

/or some universal constants M 2 > 1, ^3 G (0,1), Co, m; and d m as in 
Corollary UTt 



Proof LetC = %gP + 2. 

We proceed by induction. The case k = 1 is the result of Corollary 14.51 
with r = 1 if Co < d m , M 2 > M™ and /X3 G (/xo, !)■ Assume then it is true 
for some k and lets prove it for k + 1. It means that we are assuming that 

\{u > M 2 fe+1 } n Q x x [0, 1]| > > ^ 3 |{n > M 2 fc } n Q 1 x [0, 1]|. 

Let the dyadic disjoint covering of A = {u > M 2 fc+1 } n Qi x [0, 1] 

with respect to the fraction Lets fix m sufficiently large such that: 

(1) M2 := < I, 

(2) d m > Co, 

(3) m > rriQ with mo from Lemma 14.81 
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Lemma 14.71 and 14.81 tell us that for every piece Kj of the covering, 

\{u > M% +1 }nKj\ > m\Kj\, 

\{u > Af 2 fe+1 } n Q 1 x [0, 1]| < /x 2 | U, (Kjn 

U 3 {K 3 ) m C {u > M*} n Q 1 x [0, d m ]. 

Therefore, 

^ +1 <|{ U >M 2 fe+1 }nQ 1 x[0,l]|, 

= /i 2 (| Uj (Kj)™ n Qi x [0, 1]| + I Uj {Kj) m n Qi x [1, d m }\), 
< fi 2 (\{u > M|} n Qi x [0, 1]| + I U, (Kj) m n Qx x [l,d m ]|), 

<M2(A*3 + I u i (^i) m n<9i x [i,d m ]|). 

By fixing // 3 = (1 + /x 2 )/2 E (0, 1) we get for some constant c, 

c/4 < m c^-rnQi x [i,d m ]|. 

Therefore there is a point (a;, t) E Uj(Kj) m with t > C/U3 + 1. Hence there is 
a dyadic box Q r (xo) x [to, to + r°"r] in the covering with 

a. - 4(m + l)r <T r > C/X3, 

b. - |{u > M 2 fc+1 }nQ r (x ) x [t , *o + ^11 >Mi|Q r (a;o) x [t ,t + r a r}\, 
The first equation implies d^^T > Co for N large enough, independent of 
k, a and r. 

In order to apply Corollary 14.51 (rescaled) we need to check at least that 
(4.2) Q x x [C - 1, C ] C U*-}Q 3 <, ir (x ) x [t + ^r CT r, t + d m r CT r] 

By having chosen Co — 2 = 8 ^^ff 1 - we get that for (xq, to) G Qi x [0,1] and 

(cc, t) E Q± x [Co — 1, Co], t — to > ^y^i ° and |sc — xo| < \fn. Therefore 

Q l x [C - 1, C ] C {r < t - t < r|6(x - x )|7(3 ff - 1)}- 

Using the previous relation and that djq^r a T > Co we verify (|4,2p . It allows 
us then to conclude that \{u > M Nk }CiQ r (xo) x [t , to + r°"r]| < /ii|C- r (xo) X 
[to, to + r a r]\ which would contradict 

\{u > M 2 fc+1 }nQ r (x ) x [t ,to + r CT r]| > m\Q r {x Q ) x [t , t + r CT r]| 

if M 2 > M^. □ 

5. Regularity Results 

The purpose of this section is to prove that solutions of (II. ip are regular. 

Holder regularity follows by proving a geometric decay of the oscilla- 
tion of the solution. Holder regularity for the spatial gradient holds if the 
equation is translation invariant and the operator is elliptic with respect to 

A. 
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Holder regularity for the time derivative does not hold even for the frac- 
tional heat equation in a bounded domain if the boundary data is not suffi- 
ciently nice. Consider u, the solution of u t - + (— A) a u = in B\ x (—1,0], 
with initial value u(-, —1) =0. The boundary data outside B\ is equal to 
u = (c(t + 1/2) + XB a \B 2 ( x ))X[-i,2fi](t) where the constant c > is chosen 
small enough so that u is a subsolution to the fractional heat equation in 
Si x (—1, 0]. By the comparison principle, we have that u > u in B\ x (— 1, 0]. 
Also u = in B\ x [—1,-1/2] by uniqueness. This shows that the time de- 
rivative of u have a jump at t = —1/2. This problem will be treated in a 
future work. 

Theorem 5.1 (Holder regularity). Let oo G (0,2) and a G (do, 2). Let 
u G C(£l x (—1,0]) such that it satisfies the following two inequalities in the 
viscosity sense with Co > 0, 

Ut ~ M C (a) U - ~ C in B 2 X ("ML 
U * - M Co(a) U - C in B 2 X (-M]. 

Then there is some a G (0, 1) and C > 0, depending only on n, A and o~q, 
such that for every (y,s),(x,t) G B l / 2 x (—1/2,0] 

\u(y,s) -u(x,t)\ ^^f n ., \ 

— -— — < C ^|«||Lac(fl aX [-l,0]) + II«I|C(-1,0;L1( W )) + <-0j 

Proof. Assume without loss of generality that (x, t) = (0, 0) and also that 
u G [—1/2,1/2] in B 2 x [— 1,0]. A truncation of u also satisfies a similar 
equation with a controlled right hand side. Let v = uxb 2 > then 

v t ~ M c a v > "Co - C||u|| L i (w ) in Bi x (-1,0], 
v t~ M^v <Cq + C\\u\\ l i {ui) in B x x (-1,0]. 

Assume without loss of generality that Co + C^||it||c(-i ) 0;£ 1 M) < £ o for some 
Eq > small enough to be fixed. The idea of the proof is to construct an 
increasing sequence and a decreasing sequence Mf,, such that M^ — m^ = 
2(l/2) afc and t> is trapped between m k and M fc in Bfy 2 \k x [-(l/2)°" fc , 0]. 

Take initially m,Q = — 1/2 and Mo = 1/2. Assume we have constructed 
the sequences up to some index k. We want to find now how to construct 
77ifc + i and Mfc_|_i. Consider 

_ ?;((l/2) fc x, (l/2)^ fc t) - (m fc + M fc )/2 

so that u> G [—1,1] in Bi x [—1,0] and satisfies 

wt-^u; > -e 2- (,J - Q)fc > -e in ^ x (-1,0], 
«7t - A^ Q w < e 2- {a ~ a)k < e in £ a x (-1,0], 

given that a G (0, ctq)- 
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In B1/2 x [—1/2,0], w is between —1 and 1 and either is above or below 
zero at least 1/2 of the measure of Bin x [—1,-1/2]. Assume, without loss 
of generality, that 

|{^>0}nff 1/2 x [-i,-i/2]| > i 

1^1/2 x [-1,-1/2] I "2' 

Then we want to use Theorem 14.11 to show that w + 1 raises from in 
B\j2 x [—1/2, 0]. Still w can have negative values outside B\ which will not 
allow to apply such theorem. However, w + 1 > — 2(|2x| a — 1) outside B\, 
so that w + satisfies in -B3/4 x (—1,0], 

(w + 1)+ - M Co (w + 1) > -e - ||2(|2x| a - 1)xbl»\ Bi \\ lHu) , 

By choosing a small enough we can get right hand sides of magnitude 2eq. 
By Theorem 14. II we get that, 

n( ■ f ^9 ^ >> l^i/2 x hi/2, 0]| 
C mf w + 2e > , 

\B 1/2 x [-1/2,0] J 2 

Which for eo small enough implies that inf B 1/2 x [-1/2,0] w > 9 for some uni- 
versal 9 > 0. 

Then in this case we can choose Mk+i = and mt+i = + 9{Mk — 
m,fc)/2. Notice that, M^+i — mk+i = (1 — 9/2)2~ ak which can be made 
smaller than 2~ ak for a even smaller. In the case that 
|{^>0}n£? 1/2 x [— 1, — 1/2]| > 1 
\Bi/2 x [-1,-1/2] I -2' 

we arrive to the conclusion that M^+i = Mk — 9(Mk—mk)/2 and rrik+i = rrik 
satisfies all the inductive hypothesis which concludes the proof. □ 

Theorem 5.2 (C 1,a Regularity for translation invariant operators). Let 
(To G (0,2), a G (o"o, 2) and / G C([— 1,0]). There is p depending only on 
n, A, A, Co, suc/i i/iai ?// is an inf sup (or sup inf) type operator, translation 
invariant in space and elliptic with respect to C\{a, K, po) and u G C{B\ x 
[—1,0]) is a viscosity solution of the equation, 

Ut — Iu = f{t) in B2 x ( — 1, 0], 

then u is C 1,a in space for some universal a G (0, 1). More precisely, there 
is a constant C > 0, depending only on n, A and a, such that for every 
(x,t),(y,s) G B 1/4 x (-1,0] 

\u Xi (x,t)-u Xi (y,s)\ . n f l , 11 

(\ x - y \ + \ t - s \l/°)<* ~ {" U "L°°(B 2 x[-lfi]) + ll«IIC(-l,0;£i(w)) 
+ II/IIl°°((-i,0])) , for i = 1,. . . ,n. 
Proof. Assume without loss of generality that 

IMIz,°°(B 2 x[-:l,0]) + II/IU°°((-1,0]) < 1- 
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By considering the cutoff of u in Bi we have, as in the previous proof, that 
we can assume also that u = in W 1 \ B2. 

Let a the Holder exponent obtained by Theorem 15.11 and assume that 
it is not the reciprocal of an integer by making it smaller if necessary. Let 
5 = l/(4|l/aj). Fix a unit vector e G W 1 , a number h E (0,5/8) and, 
for k = 1, 2, . . . , [l/a\ , let rj k be a smooth cut-off function supported in 
B(3/4-k8)-8/4 an d equal to one in Bt 3 u_ k g\_g/ 2 . Define the following incre- 
mental quotients, 

hk/ . u(x + he, t) — u(x, t) 
w n ' k (x,t) = — 



ak 



h,k, ,s (r) k u){x + he,t) - (?] k u)(x,t) 



\h\ 



ak 



h,k, ((1 - V k )u)(x + he,t) - ((1 - y k )u)(x,t) 
w 2 = \^ETk • 

We will show that for every k = 0, 1, . . . , |_X /ckJ — 1, if 

\\ wk,k \\L°°(B 3/i _ kS x{-(3/4-kS),0]) < C{k) 

then in B 3 / 4 _( k+ i^s x [—(3/4 — (k + 1)5), 0] the following estimate holds 

(5 - 2) (| X _ y | + | t - S |l/^ <^(*+lJ 



where w h, ° = u and the constants C(k) are independent of h. 

When (x,t) £ B( 3 /4- fe< 5)-5/8 x [-((3/4 - kS) - 6/8), 0], |^' fc | is bounded 
above by the product rule and the hypothesis (|5.1|) . 

\w^ k (x,t)\ < C{k) + ||»7 fc |U«'(B ( a /4 _ H) _ 4/8 x[-((3/4-fc«)-*/8) J 0])- 

If x e (M n \ S (3/4 _ fc<5) _ 5/8 ) x [-((3/4 - M) - 5/8), 0] then ^(s.f) just 
cancels. 

By using that the equation is translation invariant we have that u and 
u(- + he, ■) satisfy equations in the same ellipticity family. By Theorem 12.21 
w h,k satisfy the following inequalities in B\ x (—1,0] in the viscosity sense, 

w^ k -M- c w h ' h >0, 
w^ k -M^w h ' k <0. 

The function satisfy a similar equation as w h,k in -B(3/4_fc<5)_3<5/4 x 
(—1,0], the difference is on the right hand side introduced by the cutoff, 



/ h,k\ » A ~ h,k ^ 1 A — h,k 

{w{ )t-M c w{ >-M c w 2 , 

1 h,k\ 1 A + h,k ^ 1 v+ h,k 
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For x G -B(3/4-fc<5)-3<5/4 the terms are controlled by 

by using that 

\K[ y )-K(y-h)\ 



j 



\h\ 



-dy < C every time \h\ < 



Po 



with po = S/8. Indeed, for L G C\ with kernel K, x £ £?(3/4-fc5)-3,5/4, 



y\ < 5/8 we have that W2' k (x + y) = and by the product rule 



Lw^ k {x)\ 



W2 h (x + y,t)K(y)dy 

(1 — rj k )u(x + y + h) — (1 — r] k )u(x + y) 



(1 -?7 fc )u(x + y)|/ l | 



l^afc 

^akKM-Kiy-h) 



\h\ 



K{y)dy 
dy 



c C. 



We get then the equations for w^' in -6(3/4- ^5) -35/4 X (—((3/4 — kS) — 
35/4), 0] 

(^ fc ) 4 -^^< fc >-C, 
(^ fc ) t ->[+^ fc <C7. 

By applying Theorem 15.11 to from B^/i-kS)-^/^ x ( — ((3/4 — kS) — 
35/4), 0] to i?3/4_(fc + i)5 x (— (3/4— (/c+l)<5), 0] we conclude that for a constant 
C{k + 1) independent of h, 



\Wi (M) - w{ (y,s) 



< C(k + 1). 



(\x -y\ + \t- s\ l l a ) a 

Which is equivalent to (|5.2p . 

By Lemma 5.3 in [I] we get that (|5.2p implies that w h,k+1 is also bounded 
by a constant independent of h. Therefore we can apply this procedure up 
to obtaining that u x is bounded in -B 1 / 2 x [—1/2,0]. Finally, by applying 
the previous argument one more time to the Lipschitz quotient we conclude 
the theorem. □ 
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